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1.1 Introduction 

Thermal fluctuations alone cannot create a steady directed transport in an 
unbiased system. However, if a system is out of equilibrium, the Second Law of 
Thermodynamics no longer applies, and then there are no thermodynamical 
constraints on the appearance of a steady transport [TJ [5] . A directed current 
can be generated out of a fluctuating (time-dependent) external field with zero 
mean. The corresponding ratchet effect [21 HI El [8] has been proposed as 
a physical mechanism of a microbiological motility more then a decade ago 
[11 [5] . Later on the ratchet idea has found diverse applications in different areas 
[51 [71 [H] , from a mechanical engine [9] up to quantum systems and quantum 

devices nniinj [m[T3J El ei . 

When the deviation from an equilibrium regime is small (the case of weak 
external fields) one may use the linear response theory in order to estimate 
the answer of the system p~6l Ej. However, due to the linearization of the 
response, the current value will be strictly zero since the driving field has 
zero bias. Therefore, one has to take into account nonlinear corrections and 
then derive the corresponding nonlinear response functional |17[ 1 1 8] , which 
may become a very complicated task, if the nonadiabatic regime is to be 
considered. 

To obtain a dc-current, one has to break certain discrete symmetries, which 
involve simultaneous transformations in space and time. A recently elaborated 
symmetry approach [19j j20] established a clear relationship between the ap- 
pearance of a directed current and broken space-time symmetries of the equa- 
tions of motion. Thus, the symmetry analysis provides an information about 
the conditions for a directed transport appearance without the necessity of 
considering a nonlinear response functional. 

Most theoretical and experimental studies have focused on ratchet real- 
izations at a noisy overdamped limit [SEllE]. However, systematic studies of 
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the underlying broken symmetries, and the largest possible values of directed 
currents achieved for different dissipation strength, show that the dc current 
values typically become orders of magnitude larger in the limit of weak dissipa- 
tion [211 122] . The corresponding dynamics is characterized by long space-time 
correlations which may drastically increase the rectification efficiency [22] [23] . 

Fast progress in experimental studies of cold atoms ensemble dynamics 
have provided clean and versatile experimental evidence of a ratchet mecha- 
nism in the regime of weak or even vanishing dissipation [24, 25]. The results 
of the corresponding symmetry analysis for the regime of classical dynamics 
has already been successfully tested with cold Rubidium and Cesium atoms 
in optical lattices with a tunable weak dissipation [26] [27] . Further decreasing 
of the dissipation strength leads to the quantum regime [24] . Recent exper- 
iments have shown the possibility to achieve an optical lattice with tunable 
asymmetry in the quantum regime 28J. A Bose-Einstein condensate (BEC) 
loaded into an optical potential is another candidate for a realization of quan- 
tum ratchets in the presence of atom-atom interactions [25]. While there is 
obvious interest in experimental realizations of theoretically predicted sym- 
metry broken states, another important aspect of the interface between cold 
atoms and the ratchet mechanism is, that new possibilities for a control of the 
dynamics of atomic systems by laser fields may be explored GU] ■ 

The objective of this work is to provide a general introduction into the 
symmetry analysis of the rachet effect using a simple, non-interacting one- 
particle dynamics. Despite its simplicity, this model contains all the basic 
aspects of classical and quantum ratchet dynamics, and may be used also as 
a starting point of incorporating atom- atom interactions. 

1.2 Single particle dynamics 

We start with the simple model of an underdamped particle with mass m, 
moving in a space-periodic potential U (x) = U(x + A) under the influence of 
the external force x(t) with zero mean: 

mx + 7 z - f(x) - x (t) = 0. (1.1) 

Here f{x) = —U'(x), J Q A f(x)dx = 0, and 7 is the friction coefficient. Next, we 
ask whether a directed transport with nonzero mean velocity, (x) ^ 0, may 
appear in the system (jl.l| . 

If x(t) = €{t) is a realization of a Gaussian white (i.e. delta-correlated) 
noise, obeying via its correlation properties the (second) fluctuation-dissipation 
theorem [17] , Eq. (|1.1[) then describes the thermal equilibrium state of a parti- 
cle interacting with a heat bath. From the Second Law of Thermodynamics it 
is follows that a directed transport is absent, independently of the particular 
choice of the periodic potential U(x) [4] 16] 121], 
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The presence of temporary correlations in x(t) may change the situation 
drastically. A simple way to get such correlations is to use an additive periodic 
field E(t), 

X(t)=m+E(t) , E(t)=E(t + T) , f E(t)dt = 0. (1.2) 

Jo 

If £(t) is a realization of a white noise, then the functions — and 
£(t + r) are also realizations of the same white noise, and their statistical 
weights are equal to the statistical weight of the original realization. For what 
comes, the noise term £(£) will thus not be relevant for the following symmetry 
analysis. We consider the symmetries of the deterministic differential equation 

mx + jx-f{x) -E{t) =0. (1.3) 

Eq. (|1.3p contains two periodic functions, f(x) and E(t), both with zero mean. 
The properties of the symmetries of the Eq. l|1.3[) are strongly depending on 
the symmetry properties of these functions. 

1.3 Symmetries 

1.3.1 Symmetries of a periodic function with zero mean 

Let us consider a periodic function g(z + 2n) = g(z) with zero mean, 
/ T g(z)dz = 0. This function can be expanded into a Fourier series 

oo 

g(z) = ^2 9k-exp(ikz), (1.4) 

k— — oo 

where go = 0. We will consider real- valued functions; therefore, gu = .gl^. 

The function g{z) may possess three different symmetries. First, it can 
be symmetric, g(z + zq) = g(—z + zq), around a certain argument value zq. 
For such functions we will use the notation g s . The Fourier expansion (|1.4[) 
contains, after the shift by zq, only cosine terms, so gk(zo) = 9k ■ exp(ikzo) 
are real numbers, i.e. gk(zo) — g-k(zo)- 

Second, the function g(z) can be antisymmetric, g(z + zi) — — g{— z + z\), 
around a certain value of the argument, z±. For such functions we will use 
the notation g a . The corresponding Fourier expansion (jl.4p contains only sine 
terms (after the shift by z\), so gk(zo) — gk • exp(ikzo) are pure imaginary 
numbers, and gk(zo) = -g-k(zo). 

Finally, the function g(z) can be shift- symmetric, g(z) = —g(z + it). The 
Fourier expansion of a shift-symmetric function g s h{z) contains odd harmonics 
only, g 2m = 0. 

It is straightforward to show that a periodic function g(z) can have cither 
none of the above mentioned symmetries, or exactly one of them, or all three of 
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them. Let us consider several simple examples. The function cos(z) possesses 
all three symmetries. The function cos(z) + cos(3z + <f)) always possesses shift- 
symmetry and in addition may be simultaneously symmetric and antisymmet- 
ric for = 0, ±7r. The function cos(z)+cos(2z + 4>) is not shift-symmetric, thus 
it will either have no other symmetry at all, except for = 0, ±7r (symmetric), 
and 4> = ±7r/2 (antisymmetric). 

1.3.2 Symmetries of the equations of motion 

The system dynamics in eq. (jl.3p can be described by three first-order au- 
tonomous differential equations, 

x = - , P = f{x)+E(T)-2-p, f = l. (1.5) 
m m 

The phase-space dimension is three. We are looking for symmetry transfor- 
mations, S, which do not change the equation (|1.5|l . but do change the sign 
of the velocity x. Such transformations map the phase space {x,p, r} onto 
itself. If we find such a transformation, we then apply it to all points of a 
given trajectory. We get a new manifold in phase space, which also represents 
a trajectory, i. e. a solution of the equations (jl.5p . The original trajectory and 
its image may coincide (or may not). 

Let us assume that we have found such a transformation. Next, we consider 
the mean velocity, v = lim^oo ^(x(t + s) — x(t )), on the original trajectory. 
If the trajectory and its image coincide, then v = 0. If they are different 
then their velocities have the same absolute value but opposite signs. If, in 
addition, both the trajectories have the same statistical weights in the presence 
of a white noise, then we can conclude that the average current in the system 
(]1.3|) is equal to zero [19]. 

There are only two possible types of transformations which change the 
sign of the velocity x: they include either a time-reversal operation, t — > — t, 
or a space inversion, x — > —x (but not both operations simultaneously!). 

The following symmetries can be identified |19j : 

T 

S a ■ x -x , t-^t + — , if {f a , E sh } ; 

S b : x -» x , t -» -t , if {E s , 7 = 0} ; (1.6) 

S c ■ x-*x + -,t->-t, if {f s h , E a , m = 0} . 

The symmetry Sb requires zero dissipation, 7 = 0, i.e. it requires the Hamil- 
tonian regime, and the symmetry S c can be fulfilled in the overdamped limit 
(i.e. m = 0) only. Note that all symmetries require certain symmetry proper- 
ties of the function E(t). Usually, an experimental setup allows to tune the 
shape of the time-dependent field E(t) easier than the shape of the spatially 
periodic potential [21] • A proper choice of the force E(t) may break all three 
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symmetries for any coordinate dependence of the force f(x). We restrict the 
further consideration to the case of a symmetric potential U (x) = 1 — cos(a;) 
while using a bi-harmonic driving force, 

E(t) = E 1 cos{t)+E 2 cos(2t + 6), (1.7) 

for a symmetry violation. If 9 ^ 0, it/2, tt, 3it/2 then all three symmetries (jl.6[) 
are broken and we may count on a nonzero mean velocity, v ^ 0. 

1.3.3 The case of quasiperiodic functions 

We generalize the symmetry analysis to the case of quasiperiodic driving field 
E{t) [31 133] • 

We consider a quasiperiodic function g(z) to be of the form 

z ' 

g{z) ee g(zi, z 2 , zn) , -^j = &i (1-8) 

where all ratios flij ' flj are irrational if i ^ j and g(zx, z 2 , z% + 27r, zn) = 
g(z%, z 2 , Zi, zpf) for any i. Such a function may have numerous symme- 
tries. With respect to the following symmetry analysis of the equation of 
motion we will list here only those symmetries of g which are of relevance. 
It can be symmetric g s {z\ 1 z 2 , z^) = g s (—zi, —z 2 , — zn), antisymmetric 
(ja(zi, z 2 , zpf) — —g a (—Zt, —z 2 , — zn). It can be also shift-symmetric for 
a given set of indices g s h,{i,j,—,m} which means that g changes sign when a 
shift by 7r is performed in the direction of each zi, Zj, z m only, leaving the 
other variables unchanged. 

The relevant symmetry properties of g are thus studied on the compact 
space of variables {zi, z%, Zn}- The irrationality of the frequency ratios 
guarantees that in the course of evolution of z this compact space is densely 
scanned by these variables with uniform density in the limit of large z. At the 
same time we note that it is always possible to find a large enough value Z 
such that 

1 r 

lim - / (g(z + Z) - g{z)f&z < e (1.9) 

T^OO T J Q 

with (arbitrarily) small absolute value of e. For a given value of e this defines 
a quasiperiod Z of the function g(z). 

In order to make the symmetry analysis of the equation of motion trans- 
parent, we rewrite it (skipping the noise term) in the following form |33j : 

mi + 7± - f(x) - E(<f>i,<f>2, (j>N) = , (1-10) 

(pi = wi , 

4>2 — U 2 , 



4>N = U N 
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The function f(x) is also assumed to be quasiperiodic with M corresponding 
spatial frequencies. 

The following symmetries can be identified, which change the sign of (x) 
and leave (jl.lOp unchanged: 

■ X ^ X : 4*i,j,...,m ^ 0jj,...,m ~t~ ^ ? if {fa i -^s/i.{z,j, . . . ,m} } 3 

S b : x -> x , t -> -t , if {E s , 7 = 0} ; (1.11) 

S c : x — ► x + | , t — > -t , if {/ s h,{i,2,3,...,M} i E a , m = 0} . 

The symmetry S a is actually a set of various symmetry operations which are 
defined by the given subset of indices {i, j, ...,m}. 

The prediction then is, that if for a given set of parameters any of the 
relevant symmetries (|1.11[) is fulfilled, the average current will be zero. If 
however the choice of functions f(x) and E(t) is such that the symmetries are 
violated, a nonzero current is expected to emerge. 

1.4 Dynamical mechanisms of rectification: The 
Hamiltonian limit 

Let us consider the limit 7 = (Hamiltonian case) [THl Due to time 
and space periodicity of the system (jl.3p we can map the original three- 
dimensional phase space (x,p,t) onto a two-dimensional cylinder, T 2 = 
(a; modl,p), by using the stroboscopic Poincare section after each period 
T = 2tt/lu. For given initial conditions {x(0),p(0)}, we integrate the system 
over time T, and then plot the final point, {x(T),p(T)}, on the cylinder T 2 . 

For E{t) = the system (|1.3[) is integrable and there is a separatrix in the 
phase space which separates oscillating and running solutions. A non-zero field 
E(t) destroys the separatrix and leads to the appearance of a stochastic layer 
(see Fig. ll . 1[) . In this part of the phase space the system dynamics is ergodic, i. 
e. all average characteristics are the same for all trajectories, launched inside 
the layer. Therefore, the symmetry analysis is valid for all trajectories on 
this manifold. Numerical studies have confirmed this conclusion [T9l l23l [20] , 
Fig. II. 21 shows several trajectories x(t) from chaotic layers and illustrates the 
fact that the violation of symmetries causes a directed motion of the particle. 

The dynamics within the stochastic layer can be roughly subdivided into 
two distinct fractions. The first one corresponds to ballistic flights near the 
layer boundaries. They appear due to a sticking effect [53]. A random diffu- 
sion within a chaotic bulk is attributed to the second fraction. A rectification 
effect appears due to a violation of the balance between ballistic flights in op- 
posite directions [23] . This interpretation supports the view, that even in the 
presence of damping and noise, the ratchet mechanism relies on harvesting on 
temporal correlations of the underlying dynamical system. Ballistic flights are 
just such examples of long temporal correlations on a trajectory which is over- 
all chaotic. Therefore it is not surprising, that the ratchet effect is stronger in 
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x 



Fig. 1.1. Poincare map for the system (|1.5ll.7p . The parameters are Ei — 0.252, 
E 2 = 0.052, 7 = 0. (a) 6 = 0; (b) 6 = tt/2. 

the dissipationless limit, since dissipation will introduce finite (and possibly 
short) time scales which cut the temporal correlations down. The averaged 
drift velocity can be estimated by using a sum rule [35] . From the correspond- 
ing approach, which is based on a statistical argument by the authors of Ref. 
|35j . it follows, that a mixed space, i.e. a stochastic layer with boundaries and 
embedded regular submanifolds (islands), presents the necessary condition for 
a directed transport. 

The adding of a non-zero dissipation, 7^0, does not change the situation 
drastically [20]. The symmetry analysis is still valid for this case. The phase 
space is shared by different transporting and non-trasnporting attractors with 
their corresponding basins of attraction, which are strongly entangled inside 
the former stochastic layer region. A symmetry violation causes a desym- 
metrization of basins. Finally, a weak noise leads to a trajectory wandering 
over different basins, sticking to corresponding attractors, and, finally, to the 
rectification effect. The long flights which appear at the Hamiltonian limit 
are damped after a characteristic time which is the shorter, the larger the 
dissipation strength 7 is [23l [20] . 

A systematic analysis shows that, under the condition of full symmetry vi- 
olation, the approach of the dissipationless limit leads to a drastic increase of 
the dc current value [22] , which depends on the characteristics of the stochas- 
tic layer [23j . It has been shown that, in a full accordance with the symmetry 
analysis, the dc current disappears near 9 = 0, it for the case of weak dissi- 
pation, and near 9 = ±/2 at the strong dissipation limit. The value of the 
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phase 0, at which the current becomes zero, is a monotonous function of the 
dissipation strength 7 [22] . 

An inclusion of a dc-component to the external field, E(t) = E(t) + E^c, 
may lead to a directed transport against a constant bias E^ c , even in the 
Hamiltonian limit [36] . 

The abovementioned results have been confirmed in cold atoms experi- 
ments, performed in the group of Renzoni [26] , In these experiments, atoms 
of Cs and Rb have been cooled to temperatures of several mK. An optical 
standing wave, created by a pair of counter-propagating laser beams, formed 
a periodic potential for the atoms. Finally, a time-dependent force E(t) has 
been introduced through a periodic modulation of the phase for one of the 
beams. The results of the above symmetry analysis have been verified by 
changing the relative phase <p and by tuning the effective dissipation strength. 

The case of the quasiperiodic driving force E(t) for cold atoms ratchets 
also has been studied experimentally [27], with a similar outcome. 



1.5 Resonant enhancement of transport with quantum 
ratchets 

A quantum extension of the (dissipationless) system dynamics in eq. (|1.3j) 
can readily be achieved [37]. The system evolution can be described by the 
Schrodinger equation, 
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ih ¥t m))=m)m)), (i.i2) 

where the Hamiltonian H is of the form 

2 

H(x,p, t) = y + (1 + cos(x)) - xE(t). (1.13) 

The system (jl . 12[) describes a cloud of noninteracting atoms, placed into a pe- 
riodic potential (formed by two counter-propagating laser beams) and exposed 
to an external ac field fTTZjU . 

Because of the time and space periodicity of the Hamiltonian (|1.13p , the 
solutions \if) a (t + to)) = U(t,to)\ip a (to)} of the Schrodinger equation (|1.12p 
can be characterized by the eigenfunctions of the Floquct operator U(T,to) 
which satisfy the Floquet theorem: \i() a {t)) — exp(— i^t)\(f) a (t)), \<j) a (t+T)) = 
\<p a (t)) (here to is the initial time). The quasienergies E a (— 7r < E a < n) and 
the Floquet eigenstates can be obtained as solutions of the eigenvalue problem 
of the Floquet operator 

U(T,t )\<fi a (t )) = e- lE -\<f> a (to)) (1.14) 

with a denoting the band index and with k being the wave vector [?. I39[[38] . 
An initial state can be expanded over Floquet-Bloch eigenstates, \ip(to) = 
Eai^a,fc(^)|^a,k) and the subsequent state's evolution is encoded in the 
coefficients {C a .k}- We restrict further consideration to the case n = which 
corresponds to initial states where atoms equally populate all (or many) wells 
of the spatial potential. 

The mean momentum expectation value, 

J(t ) = lim l/t f (ip(t,to)\p\Tp(t,t )), (1.15) 

measures the asymptotic current. Expanding the wave function over the Flo- 
quet states the current becomes 

J(to) = J2 { J 3 )"\ C M\ 2 ' (1-16) 

a 

where (p) a is the mean momentum of the Floquet state \cj> a ) [37l [38] , 

The analysis of the transport properties of the eigenstates shows that the 
quantum system inherits the symmetries of its classical counterpart [37j . In 
particular, the symmetries of the classical equations of motion translate into 
symmetries of the Floquet operator. The presence of any of these symmetries 

3 The dissipation may be included into quantum dynamics by coupling the sys- 
tem (|1.13|l to a heat bath, Hdiss{x,p, t, {q}) = H(x,p,t) + Hb{x, {q}). Here 
Hb{x, {q}) describes an ensemble of harmonic oscillators {q} at thermal equilib- 
rium interacting with the system [TO] , 
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Fig. 1.3. (Color online) (a) Poincare section for the classical limit, (|1. 7111. 13]) ; (b-f) 
Husimi representations for different Floquet eigenstates for the Hamiltonian f|l . 13f) 
with h — 0.2 (momentum is in units of the recoil momentum, p r = %kh, with 
fc_L = 1). The parameters are E\ — E2 = 2, ui — 2, 9 — —n/2 and to = for (b-e), 
and E 1 = 3.26, E 2 = 1, (J = 3, 9 = -ir/2 and to = for (f). 



results in a vanishing the time-averaged expectation value of the momentum 
operator for each Floquet eigenstate: (p) a — 0. Thus, if one of the symmetries, 
S a , Sb (jl.lip . holds then (p) a — for all a. Consequently J(t ) — in this 
case. 

By using the Husimi representation [40j we can visualize different eigen- 
states in the phase space, {x,p, t} and establish a correspondence between 
them and the mixed phase space structures for the classical limit (Fig. ll.3|) . 

Since the Schrodinger equation (|1.12[) is linear, the system maintains a 
memory of the initial condition for infinite times [41] . The asymptotic current 
value depends on the initial time, to, and on the initial wave function, ip(ta). 
For a given initial wave function, \ip) = |0), we can assign a unique current 
value by performing an averaging over the initial time to, J = l/T J Q J(to)dto 
[37] . Fig |1.4l shows the dependence of the average current on the asymme- 
try parameter 9. Sharp resonant peaks for E2 = 1.2 where the current value 
changes drastically are associated with interactions between two different Flo- 
quet eigenstates. The Husimi distributions show that one state locates in the 
chaotic layer, and another one in a transporting island. Off resonance the ini- 
tial state mainly overlaps with the chaotic state, which yields some nonzero, 
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Fig. 1.4. (a) The average current J (in units of the recoil momentum) vs 9 for 
different amplitude values of the second harmonic, E2: 0.95 (pointed line), 1 (dashed 
line) and 1.2 (solid line). Insets: relevant details of the quasienergy spectrum versus 
6 in the resonance region for £2 = 1 (top right) and E2 = 1.2 (bottom left). The 
parameters are E\ — 3.26 and ui — 3. (b) The average current J (in units of the 
recoil momentum) vs 8 for E\ — 3, E2 = 1.5 and u = 1. 

yet small, current. In resonance Floquet states mix, and thus the new eigen- 
states contain contributions both from the original chaotic state as well as 
from the regular transporting island state. The Husimi distribution of the 
mixed state is shown in Fig ll.3f f). the strong asymmetry is clearly observed. 
The regular island state has a much larger current contribution, resulting in 
a strong enhancement of the current. 

To conclude this section, we would like to emphasize the following two 
points. For both cases, i.e. the classical and the quantum one, the overall, 
total current over the whole momentum space is zero |35j . Thus, it is essential 
to have the initial state prepared localized near the line p = 0, because for 
broad initial distributions the asymptotic current tends to zero. However, if 
the dynamics is restricted to the lowest band of the periodic potential, no 
current rectification does occur l42l . 



1.6 Outlook 

This surveyed symmetry analysis, originally put forward in Refs. [19 s 20], pro- 
vides a general toolbox for the prediction of dynamical regimes for which one 
can (or cannot) obtain the rectification and directed current phenomenon for 
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a given transport dynamics. First, one has to set up the equations of motions 
and define a the observable (current, magnetization, angular velocity, energy 
flux, etc) which should become nonzero, in terms of these dynamical variables. 
Then, one examines whether there exist transformations (symmetries) which 
change the sign of the observable and at the same time leave the equations 
of motion invariant. Upon breaking all the symmetries one can expect the 
emergence of a non-zero, directed current. This strategy has been successfully 
tested with Josephson junctions (fluxon directed motion) [43] [44] and as well 
with paramagnetic resonance experiments (spin magnetization by a zero-mean 
field) [45l|46]. 

Herein, we focused only on the one-dimensional case. By use of more laser 
beams, experimentalists can fabricate two- and three-dimensional optical po- 
tentials [25]. By changing the relative phase between lattice beams, 2D- and 
3D-potentials with different symmetries and topologies can be achieved |47j . 
This fact incites for an extension of the present ratchet studies into higher 
dimensions. 

Moreover, for the phenomenon of Bose-Einstein-condensation (BEC) of 
cold gases, interactions between atoms become essential and nonlinearities 
start to play an important role [25]. Many features of BEC dynamics are 
manifestations of general concepts of nonlinear physics, such as soliton cre- 
ation and propagation. These collective excitations can then themselves be 
subjected to a ratchet transport mechanism [48] . 
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